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Upward-going muons observed by the Super-Kamiokande detector are produced by high-energy 
atmospheric neutrinos which interact in rock around the detector. Those which pass completely 
through the detector have a mean parent neutrino energy of ~ 100 GeV, while those which range 
out inside the detector come from neutrinos of mean energy ~ 10 GeV. The neutrino baseline varies 
with the observed muon zenith angle, allowing for an independent test via disappearance of the 
neutrino oscillations observed in the Super-Kamiokande contained events. 614 upward through- 
going and 137 upward stopping muons were observed over 537 (516) live days, resulting in a flux of 
$t = 1.74±0.07(stat.)±0.02(sys.) ($ s = 0.380±0.038(stat.)±°'°^(sys.)) xl0 _13 cm- 2 s _1 sr _1 . The 
observed stopping/through-going ratio 1Z = 0.218 ± 0.023(stat.)^p Qjg(syst.) is 2.9a lower than the 
expectation of 0.368lQ 044(theo.). Both the shape of the zenith angle distribution of the observed flux 
and this low ratio are inconsistent with the null oscillation hypothesis, but are compatible with the 
previously observed z/ M <-> v T oscillations. Taken as a whole, the addition of these higher energy 
data to the contained neutrino events provides a better measurement of the oscillation parameters, 
narrowing the allowed parameter range to sin 2 28 > 0.9 and 1.5 x 10 _3 eK 2 < Am 2 < 6 x 10 -3 at 
90% confidence. 

I. INTRODUCTION 

The Super-Kamiokande ( "Super-K" ) atmospheric neutrino analyses presented in |Q and updated in these proceed- 
ings H| provide evidence for <-> v T flavor oscillations. These analyses examine neutrinos whose interaction vertices 
are inside the water of the Super-K detector However, to extend this analysis to higher energies, where the 

steeply falling power law spectrum of the neutrinos' cosmic ray parents reduces the absolute flux of atmospheric 
neutrinos, a larger interaction volume must be used. This is accomplished by studying the products of neutrino 
interactions in the rock surrounding the detector. 

Energetic atmospheric or interact with the rock surrounding the Super-K detector and produce muons via 
weak interactions. While the constant rain of cosmic ray muons overwhelms any downward-going neutrino induced 
muons, upward-going muons are neutrino induced because upward-going cosmic ray muons cannot penetrate the 
whole Earth. The flavor of these parent neutrinos is or v^, as z/ e and z7 e induced electrons and positrons shower 
and die out in the rock before reaching the detector. Detection of upward-going muons resulting from from tauons 
produced in v T or v T charged current interactions is suppressed by branching ratios and kinematics to < 3% of the 
induced muon flux. Thus, the upward-going muon technique results in a reasonably pure sample of muon-flavored 
neutrinos, allowing a test of possible disappearance due to flavor oscillations. 

Those muons energetic enough to cross the entire detector are defined as "upward through-going muons" Q . The 
mean energy of their parent neutrinos is approximately 100 GeV. Those upward-going muons that range out inside 
the detector are defined as "stopping upward-going muons" , and come from parent neutrinos with a mean energy of 
about 10 GeV. In comparison, the "Partially Contained" events from Q have a similar parent neutrino spectrum to 
that of the stopping upward-going muons, the "Multi-GeV" neutrinos from typically have several GeV of energy, 
and the "Sub-GeV" neutrinos from j^| are less than a GeV. The upward-going muon parent neutrino spectra are 
shown in Fig. [l]. 

II. THE EXPERIMENT 

The Super-K detector is a 50 kton cylindrical water Cherenkov calorimeter located at the Kamioka Observatory 
and administered by the Institute for Cosmic Ray Research of the University of Tokyo. To reduce the cosmic ray 
muon background, the detector is placed ~1000 m underground in the Kamioka mine, Gifu prefecture, Japan. The 
detector is divided by an optical barrier instrumented with photomultiplier tubes ("PMT"s) into a cylindrical primary 
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detector region (the Inner Detector, or "ID") and a surrounding shell of water (the Outer Detector, or "OD") which 
allows the tagging of entering and exiting particles. Details of the detector can be found in Ref. ||. An upward-going 
muon is defined as an event that appears to enter from the rock and is traveling upwards. Thus, PMT activity in the 
OD at the muon's entrance point is required, and those events which also have an OD signal at the muon's exit point 
are classified as through-going. Note that a neutrino interaction inside the water of the OD itself also produces an 
"entering" signal, so some small fraction (< 1.5%) of the upward-going muon sample actually originates in the OD 
rather than the rock. This effect is accounted for in the expected flux calculations. 

Parent Neutrino Energy Distribution (muon track length > 7m) 
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FIG. 1. The energy spectra of the parent neutrinos of stopping upward-going muons (left) and upward through-going muons 
(right). The dashed lines are the result of using the Bartol input fluxes, and the solid lines are from those of Honda. 

The cosmic ray muon rate at Super-K is 2.2 Hz. The trigger efficiency for a muon entering the detector with 
momentum more than 200 MeV/c is ~100% for all zenith angles. The nominal detector effective area for upward- 
going muons with a track length > 7m in the ID is ^1200 m 2 . 

The data used in this analysis were taken from Apr. 1996 to Jan. 1998, corresponding to 537 days of detector 
livetime for the through-going muon analysis and 516 live-days of the stopping muon analysis. Event reconstruction 
is made by means of the charge and timing information recorded by each hit PMT. The direction of a muon track is 
first reconstructed by several automated grid search methods, which find the track by minimizing the width of the 
residual distribution of the photon time-of-flight subtracted ID PMT times. Finally, an independent double hand-scan 
of all upward-going muon candidates (event loss probability < 0.01%) is then done to eliminate bad fits and to obtain 
a final precision fit. 

A minimum track length cut of 7m (~1.6 GeV) was applied. This cut serves to eliminate short pathlength events 
that are very close to the PMTs and thus hard to reconstruct, as well as providing a large energy threshold and 
eliminating non-muon showering background. To reduce the abundant downward-going cosmic ray muons, events 
satisfying cos© < 0.1 are selected, where O is the zenith angle of the muon track, with cosO < corresponding 
to upward-going events, and the down-going muons satisfying < cosO < 0.1 are used to estimate the up/down 
separation resolution near the horizon. After a visual scan by two independent groups and a final direction hand-fit, 
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614 through-going and 137 stopping upward-going muon events with cos < remain. 

Due to the finite angular resolution (1.5°) and multiple Coulomb scattering in the nearby rock, some down-going 
cosmic ray muons may appear to have cos < 0. The estimation of this background in the nearly horizontal zenith 
angle bin is done by fitting the nearly horizontal down-going cosmic ray muon zenith angle shape and projecting its 
tail below the horizon (see || for a detailed discussion). This background is estimated to be 4.3±0.4 through-going 
and 13.2±3.5 stopping events, all contained in the last (closest to horizontal) zenith angle bin. The stopping muon 
sample has a larger contamination than the through-going sample due to the lower energies allowed in the stopping 
muon sample, and is more uncertain due to lower statistics. The contamination at the Kamioka site due to cosmic 
ray photoproduced upward-going pions |6| and electromagnetic showers (from v e charged current or neutral current 
interactions) meeting the 7m track length requirement is estimated to be < 1%. 

The total detection efficiency of the complete data reduction process for upward through-going muons is estimated 
by a Monte Carlo simulation to be >99% which is uniform for — 1 < cos < 0. Using the upward/downward symmetry 
of the detector configuration, the validity of this Monte Carlo program has been confirmed by downward-going cosmic 
ray muon data. 

The resulting total flux observed in each of the two upward-going muon samples integrated over the whole lower 
hemisphere is listed in Table |. The flux plotted as a function of zenith angle can be seen in Fig. |[ The lower statistics 
available for the stopping upward-going muon analysis is the reason for the coarser binning of those events. 
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FIG. 2. Stopping (a) and through-going (b) upward-going muon fluxes observed in Super-K as a function of zenith angle. 
The error bars indicate uncorrelated experimental systematic plus statistical errors added in quadrature. The solid histogram 
shows the expected fluxes based on the Honda/GRV94 model for the null neutrino oscillation case, normalized by a = —10% 
from the best-fit no-oscillation through-going muon case. The dashed lines are the expected flux in the presence of f M «-» v T 
oscillations at maximal mixing and Am 2 = 3.2 x 10 _3 eV 2 , the best-fitting oscillation hypothesis. 



III. EXPECTED FLUX 



To predict the expected upward-going muon flux, this analysis used a model which is a combination of the atmo- 
spheric neutrino flux from Honda et al |t]] and a neutrino interaction model composed of quasi-elastic scattering || 
+ single-pion production || + deep inelastic scattering (DIS) multi-pion production. The DIS cross-section is based 
on the parton distribution functions (PDF) of GRV94DIS |l(| with the kinematic constraint of W > 1.4 GeV/c 2 . 
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Lohmann's muon energy loss formula in standard rock |Tl|] is then employed to analytically calculate the expected 
muon flux at the detector. This flux is compared to three other analytic calculations to estimate the model-dependent 
uncertainties of the expected muon flux. The other flux calculations use the various pairs of the Honda flux, the 
GRV94DIS PDF, the atmospheric neutrino flux model calculated at Bartol pj, and the CTEQ3M (k| PDF. These 
comparisons yield the theoretical systematic errors used in the calculations. In this analysis, Honda's input flux is 
used as the primary input so as to be more directly comparable with the contained event analysis JjJ, which is also 
based upon the Honda flux. 

The Honda+GRV94DIS calculation results in an expected muon flux shown in Table | for cos 9 < 0. The dominant 
error in the absolute flux comes from the absolute normalization uncertainty in the input neutrino flux, which is 
estimated to be approximately ±20% 12 7j,[l4| for neutrino energies above several GeV. The dominant error in the 
ratio comes from the uncertainty in the slope of the cosmic ray spectrum. The parent neutrino energy spectrum is 
shown in Fig. |l|. Note the clear separation in energies of the stopping and through-going upward-going muon parents. 



TABLE I. Observed and expected (Honda-GRV94) upward-going muon fluxes, in units 
of KT^cm^s^sr" 1 , and their ratio TZ = $ 3 /<l> t . 





N obs 


N exp 


Observed <3> 


Expected $ 


Through 


614 


695 


1.74 ± 0.07(stat.) ± 0.02(sys.) 


1.84 ± 0.40(theo.) 


Stop 


137 


244 


0.380 ± 0.038(stat.)+[J;°^(sys.) 


0.676 ±0.149(theo.) 


n 


0.218 ± 0.023(stat.)+^3(syst.) 


0.368^,:oS(theo.) 



IV. UPWARD-GOING MUONS AND NEUTRINO OSCILLATIONS 

The probability of a neutrino produced as a remaining that way is : 

P{v, -„„) = !- sin 2 (2tf) sin 2 ( ^f^ ) , (1) 

where Am 2 is the square of the mass difference of the two neutrino flavors, 9 is the mixing angle, L is the distance the 
neutrino has traveled since production (the "baseline"), and E v is the neutrino energy. Of these, L and E v are the 
experimental handles on the problem. While the v^'s indirectly detected via their daughter muons cannot be fully 
kinematically reconstructed in the same way as the contained neutrino interactions, there are several ways in which 
the upward-going muon signal is changed by oscillations. 

First, «-> v T oscillations produce a lower absolute flux of v^. The power of this test is reduced by the large 
(~ 20%) theoretical uncertainty in the absolute flux, but maximal mixing over many oscillation lengths would result 
in a 50% reduction in flux, so this test is still useful. 

Second, as the zenith angle of the neutrino's arrival direction changes, the baseline L changes. Note that only zenith 
angles from the nadir (cos 9 = —1) to the horizon (cos 9 = 0) are accessible in the upward-going muon analyses, as 
there is no way to tell if a downward-going muon's parent is a neutrino or, far more likely, a cosmic ray-induced muon. 
Neutrinos arriving vertically upward travel roughly 13,000 km, while those coming from near the horizon originate 
only ~500 km away. The effect of oscillations is more pronounced for neutrinos traveling the longer distances (those 
having a larger L in Eq. |l|) , distorting the shape of the zenith angle distribution in favor of fewer vertically upward- 
going muons. This shape comparison is fairly free of theoretical and systematic errors, the bin-to-bin errors ranging 
from ±(0.3 — 3.8)% ||. As can be seen in Fig. ||b, the no-oscillation through-going muon predicted shape does not fit 
the data well (x 2 /dof = 18/9), even after allowing the overall normalization to float downwards by 10%. 

Finally, since the two available samples of upward-going muons have parents of different energies, those having the 
smaller E v in Eq. Ill would have a larger chance to oscillate. This reduces the observed flux of stopping upward-going 
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muons compared to that of the higher energy upward through-going muons. Forming a ratio of the fluxes 7Z — $ s /$ t 
is a useful way to make this comparison, since much of the absolute theoretical uncertainty divides out. Thus, «-> v T 
oscillations reduce $ s , lowering 1Z. As can be seen in Table Q, the observed 1Z is 2.9a lower than that expected for 
the no-oscillations case, and Fig. ||a dramatically shows the loss of stopping muon flux. 

We estimated the most likely values of sin 2 29 and Am 2 using using both the 10 zenith angle bins of upward through- 
going muon flux and the 5 stopping upward-going muon zenith angle flux bins. The expected flux (d^>/dfl) osc for a 
given set of Am 2 and sin 2 29 is calculated and the same binning is applied to this flux as to the data. To test the 
validity of a given oscillation hypothesis, we minimize a \ 2 which is defined as: 
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where a s t a t,i (& S ys,i) is the statistical (experimental systematic) error in the observed flux (d^>/dD,Y obs for the ith 
bin, (1 + a^) is an absolute normalization factor of the expected flux, and (1 + r)) is a comparative normalization 
between stopping and through-going fluxes. The absolute flux normalization error cr Qjj is estimated to be ±22 %. a v 
is estimated to be j^j %j an d is analogous to the systematic error in the ratio 71. a sys ,i ranges from ±(0.3 — 3.8)%. 
Then, the minimum x 2 (Xmin) is sought on the Am 2 — sin 2 29 plane. 

The no-oscillation case results in a x 2 /dof = 41/15, a poor probability (3.2 x 10~ 4 ) of the null hypothesis. However, 
the best fit point of maximal mixing and Am 2 = 3.2 x 10 _3 eV 2 matches the data well, as can be seen in Fig. ||. 
Fig. [|a shows the confidence intervals in oscillation parameter space, taking into account that the overall best fit is 
slightly outside the physical region at sin 2 29 = 1.1, Am 2 = 3.7 x 10 _3 eV 2 . These results are consistent with those 
obtained in the analysis of Super-K's contained atmospheric neutrino events |Jl|pJ. However, they are obtained from 
an independent sample of much higher energy atmospheric neutrinos. 
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FIG. 3. a) the allowed region contours at 68% (dotted contour), 90% (thick solid), and 99% (dashed) C.L. obtained by 
combining the Super-K upward through-going and stopping muon data in a fit on the (sin 2 28, Am 2 ) plane for the «-» v r 
oscillation hypothesis. The star indicates the best fit point at (sin 2 28, Am 2 ) = (1.0,3.2 x 10~ 3 eV 2 ). Also shown is the 
allowed region contour (thin solid) at 90% C.L. by the Super-K contained event analysis, b) The same confidence intervals are 
presented for the global oscillation parameter fit of the Super-K contained data as well as the upward-going muons. The star 
again represents the best fit point (which happens to be at the same point in parameter space as the upward-going muon fit 
alone). The allowed regions are to the right of the contours. 
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Including both the upward-going muon terms from Eq. g and the contained event terms from and properly 
accounting for the joint systematic errors, a global oscillation fit using all of the atmospheric neutrinos available to 
Super-K can then be calculated. This global \ 2 1S of the form : 

65 5 10 5 8 , s. 2 

* 2 = E + E + E + E + EGr) ■ (3) 

FG PCfj, upthru-fi upstop-pi i=l v e4/ 

where ei/avi are the systematic error terms and their constraints. In this fit, the a flux normalization term is in 
common between the contained and upward-going muon expected fluxes, although it is not constrained, following the 
contained fit procedure. The spectral index uncertainty (E~ s from El) is common and constrained. The rest of the 
stopping/through-going uncertainty from the r\ term in Eq. |^ is placed in an 7/2 term {<j V 2 = ±7%), and a term 771 
(<7, ; i = ±7%) is introduced representing the systematic uncertainties between contained and upward-muon analyses. 
The no-oscillation case of x 2 /dof = 214.3/84 strongly rules out the null hypothesis. Minimizing the combined \ 2 
narrows the allowed region of parameter space to close to maximal mixing and 1.5 x 10 -3 < Am 2 < 6 x 10™ 3 eV 2 , 
as can be seen in Fig. ||b. This best fit point has x 2 /dof = 70.2/82. All the €{ terms fall within one sigma 
of the systematic error estimates. The overall best fit point is \ 2 /dof = 69.4/82 outside the physical region at 
sin 2 29 = 1.05, Am 2 = 3.2 x 10 _3 eV 2 , and the confidence intervals have been appropriately corrected. 

Making use of the higher energy interactions in the rock surrounding Super-Kamiokande provides an independent 
and complementary test of the «-> v T disappearance oscillations proposed to explain the atmospheric neutrino 
anomaly. These upward-going muons confirm the conclusions reached in Q, and a global fit to all the Super-K 
atmospheric neutrinos now available provides a more precise measurement of the oscillation parameters involved. 
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